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Abstract 

We propose a new algorithm, for parameter estimation that is applicable to imaging using moving and synthetic 
aperture arrays. The new method results in higher resolution and more accurate estimation than commonly used 
y—i methods when strong interfering sources are present inside and outside the field of view (terrestrial interference, 

confusing sources). 
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When multiple antennas or other sensors are used to estimate incoming signals, it is best to treat them 
all as a single array and apply one of the known array processing algorithms for the estimation. In various 
situations, however, this is not practical or even impossible. 

I n radio-astronomy, the antenna array moves with the rotation of the Earth. Correlations between the 
antenna can be obtained for a given time, but not between different points in time. Denote by R& the 
• measured correlation matrix (visibility) of the array at time t k ; the combined correlation matrix of K 
o , epochs t\, . . . , t K can be written as 
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where some of the matrix elements are unknown. The R& matrices need not be of the same size. 

The need to overcome this obstacle is as old as the days of radio astronomy itself. The standard approach 
\q to process the measured correlations (a.k.a visibility) is through using inversion (see Q]| page 128 and 
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Q where (uk,Vk) are the baselines (i.e., distances between the antennas at the measurement time), M is the 
number of measurements and 1(1, m) is the estimated power: in other words, dirty image. This is optimal 
.£h as long as epochs are independent. However, sidelobes of interfering sources are not spatially white, 
^ thus making this approach sub-optimal. The dirty image is used for further processing, by deconvolution 
algorithms to yield a better image. 

The most widely used deconvolution algorithm is the CLEAN (proposed by Hogbom [3]) and its many 
variants. The CLEAN algorithm assumes that the observed field of view is composed of point sources. 
CLEAN iteratively removes the brightest point source from the image until the residual image is noise-like. 
The point sources are accumulated during the iteration and the reconstructed image is the accumulated 
source list convolved with a reconstruction beam (usually a Gaussian). During the iterations, CLEAN 
subtracts the strongest point sources either in the image domain or in the visibility domain. The visibility 
domain CLEAN is more accurate since we are not limited to pixel resolution. Throughout this paper, we 
use the visibility domain CLEAN. Acceleration of the CLEAN algorithm can be achieved by estimating 
multiple point sources based on a single dirty image (major cycle), as well as defining windows for 
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the search procedure. Practically, defining windows reduces the size of the search space. A multi scale 
CLEAN proposed by Cornwell flU models the brightness of the sky by the sum of the components of the 
emission having different size scales. Extensions of the CLEAN algorithm to support wavelets as well as 
non co-planar arrays are reviewed by Rau et al. in [0. 

A matrix based imaging technique which proved equivalence between the radio astronomical problem 
and the signal processing formulation using beamforming was proposed by Leshem and van der Veen 
J6l and further analyzed by Ben-David and Leshem [7|. Using the matrix based technique, Leshem and 
van der Veen [8| suggested a method for the cancellation of strong interference from a small number of 
interfering sources. 



A. Matrix Based Imaging 

In this section we show the equivalence of the standard approach to dirty image calculation and classic 
(i.e., Bartlett) beamforming (see fl9), IflOl and Q for an in depth discussion). In this section we assume 
a calibrated array (extension for calibration is not hard). 

At the fc'th epoch, the measured correlation matrix is given by 
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where V(u 



is the correlation (visibility) between antenna % and antenna j located at baseline 
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on the fc'th epoch. The array steering vector is defined by 



a k (l,m) 
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where (u k , v k ) is the location of antenna i at the /c'th epoch, relative to some convenient point (u , v ), 
(I, m) are the direction cosines ,P is the number of antennas in the array and A is the wavelength. For 
K epochs, we have M = KP 2 and from Equation (|2])-((4]) we obtain that the relation between the dirty 
image and the measured correlation matrices is given by 
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which is the classic (Bartlett) beamformer with weight vector w fc (/,m) = -^a fc (Z,m). 
In the general case, for a weight vector Wfc(7, m), the dirty image is given by 



1 K 

1(1, m) = — wf (I, m)R fc Wfc(Z, m). 
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B. MVDR Beamformer 

The MVDR (Minimum Variance Distortionless Response) beam-former is designed for scenarios that 
include interfering sources in the field of view. Its weights are set to minimize the influence of the 
interfering sources while passing signals from the desired direction (i.e., to minimize the interfering 
power entering the array via its 'sidelobes'). 

The MVDR is also designed to obey the distortionless response condition; in the absence of noise, the 
array input and output signals should be equal. Thus, 

w H a = 1. (7) 
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The MVDR beamformer minimizes the total array output power. For a given observed source and 
thermal noise, the power from interfering sources is minimized. The MVDR weights are determined by 
solving the following problem: 

w = arg min w w H Rw 



(8) 



The solution to Equation ([8]) is given by 



H = — (9) 



a^R !a 

The MVDR is an adaptive method. The weights are determined by the measured visibility. This is 
unlike the classical beamformer that determines the weights according to the observing angle independent 
of other radiating sources. 



II. Adaptive Selective Sidelobe Canceller Beamforming 

In this section we present a novel image formation technique. We begin with a simple example that 
demonstrates the main idea behind the adaptive-selective-sidelobe-canceller (ASSC) algorithm; for a 
specific observation direction, the received interference through the sidelobes varies strongly as the array 
rotates. 

For simplicity, consider an East- West linear array with 20 antennas, A/2 spaced, measuring the corre- 
lations every 6 minutes for a 12-hour period. The measured correlation matrix at the fc'th epoch is R&. 
For a specific direction (l,m), the output of the fc'th beamformer, J& = w^R^w^, is composed of the 
signal-of-interest (SOI) contribution, interfering sources and the noise contribution. The contribution of the 
interfering sources is determined by their location (and strength) relative to the array sidelobes. Consider 
a scenario with a few clusters of sources (see Figure ([1^)). Figure (HJ>) shows the output power of the 
fc'th classic beamformer for a direction of a point source Si (marked at Figure <[T^)). Figure (fib) shows 
the fc'th MVDR beamformer for the same direction. From all available time epochs, only a few epochs 
yield an estimation close to the true point source intensity. The intensity estimation of most epochs is 
biased due to the interference (received through their sidelobes). The output power of the two beamformers 
(classic and MVDR) towards the empty direction S 2 (marked at Figure ([T^)), is plotted in Figures ([l]l) and 
([jj) respectively. Only a few time epochs yield a close to zero intensity estimation (the true intensity); 
whereas most of the epochs estimate biased intensity originated from the interference signal received 
through their sidelobes. The time epoch with the minimal power for a specific direction yields the best 
estimator. Averaging the output power for all epochs will result in a biased and inaccurate estimator. 

Note that although the number of reliable estimates per direction is small, the total number of correlation 
matrices for the entire FOV is much larger (it depends on the interference location and the array geometry). 
Figure ([2]) shows histograms of the number of directions (pixels) for which a specific correlation matrix 
estimated the minimal power (i.e., underwent the smallest interference). Simulation conditions are the 
same as in Figure ([T]). For the classic beamformer (Figure [2^), out of the 181 x 181 pixels in the image, 
most time epochs (more than 90%) performed best (i.e., had minimal interference) for at least 661 pixels. 
Over 65% of the time epochs benefited from minimal interference for 90% of the pixels in the image. 
As for the MVDR beamformer, of the 121 available time epochs, most time epochs (more than 90%) 
performed best (i.e., had minimal interference) for at least 507 pixels out of the 181 x 181 pixels in the 
image. Over 45% of the time epochs experienced minimal interference for 90% of the pixels in the image. 

This example demonstrates that the received array power of a specific (l,m) observation direction, 
varies significantly with the array orientation due to interfering sources. Some of the array orientations 
yield a reliable intensity estimation, whereas others yield a biased intensity estimation (aggravated by the 
interfering signals). 
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(a) Original image 
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(b) Classic beamformer outputs for Si 



(c) MVDR beamformer outputs for Si 
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(d) Classic beamformer outputs for S2 



(e) MVDR beamformer outputs for S2 



Fig. 1. The array power for Si and S2 of the classic and MVDR beamformer for all time epochs. 
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(a) Classic beamformer 



(b) MVDR beamformer 



Fig. 2. Number of pixels ((I, m) directions), a specific measurement (time epoch) estimated the minimal power (i.e., underwent the smallest 
interference). 



A. The ASSC algorithm 

Based on the observations above, we now propose an algorithm that exhibits high performance in the 
presence of interfering sources and images with a high dynamic range. For a given set of Rk correlation 
matrices, k — 1 . . . K, measured at K time epochs (or by K different arrays): 

1) Calculate the array output power (i.e., dirty image) for each epoch separately according to the desired 
beamformer, 

4(7, m) = wf (l,m)RkW k (l,m). (10) 
w fc = Wmvdr^, 7 ™) is the MVDR weight vector given by 

w^ VDR (/,m)= k \ H > k - (11) 

2) Determine the ASSC parameters k and ji k where k is the number of best epochs to consider for each 
(/,m) and fXk are their weights. These parameters are best determined using the measured data by 
plotting a histogram of the calculated Ik(l, m) for a specific (I, m). A; is determined so no epoch that 
suffers from significant sidelobes (i.e., has significantly larger power than the minimal power) will 
be selected. Typically k < 5% depending on the array geometry and the interference strength and 
location. As a rule of thumb, the stronger the interference, the smaller the A;.As for /x^, is should be 
chosen such that fik+i < /ife. 

3) For each (7, m) (each pixel in the image) find the best (i.e., smallest) k values among all measurements, 

[h(l,m),...,I~ k (l,m)] =[/ (1) , Jp),...,/^] (12) 

where [1^, 7( 2 ), . . . , L^} are the k smallest elements in the order statistics of m), . . . , Ik{1, m) 

4) Calculate the ASSC power (dirty image) according to 

k 

i ASSC (l,m) = J2^Ik(l,m). (13) 

k=l 



Similarly, the weight vector from Equation (11) can be chosen using any other beamforming technique 
(for example classic, AAR). Table (fTl) summarizes the ASSC beamformer algorithm. 
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For each incident angle (I, m): 

• Calculate the desired beamformer weight vector for each epoch. 

• Calculate the beamformer output power of each correlation matrix separately, Ik(l, m) = w^RfcWfe 

• Select the best (i.e., smallest) k measurements among m), . . . Ik(1, m) 

• Calculate the ASSC dirty image by I ASSC (l,m) = YX=i^I{k){l,m) 



TABLE I 
ASSC BEAMFORMING 




Fig. 3. Illustration of array sidelobes. 

The computational complexity of the ASSC classic/MVDR beamformer is similar Jo classic/MVDR 
beamformers respectively, with the following minor addition: for each pixel, find the k minimal powers 
from [/i(/,m) 5 . . . , I K (l, m)]. 

B. Rationale for the ASSC 

Figure [3] illustrates the sidelobes of a rotating array in two orientations, observing the same SOI (marked 
in green) in the presence of two interfering sources (marked in purple). The array at the first orientation 
(marked in red), has a strong sidelobe in the direction of the interfering sources and therefore receives 
strong interfering power. The array at the second orientation (marked in blue), receives much lower 
interference power due to the shape and location of its sidelobes relative to the interfering sources. The 
received power from the interfering sources depends strongly on the direction of the interfering sources 
relative to the array sidelobes whereas the received power from the SOI is similar for all orientations. 

The ASSC method is based on the following observations: 
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a) If a signal source is present, and noise at the antenna is neglected, all correlation matrices estimate 
the same incoming signal and its power. 

b) The different results, for different time epochs, are due to interfering signals arriving from other 
directions through the array sidelobes. 

c) By choosing the time epoch with the minimal power, we choose the correlation matrix with the smallest 
interfering power, which happens to best suppress the interference. 

Some comments are in place: 

a) The proposed method is adaptive and the selected epochs are based on the signal estimates. We 
implicitly assume that at least one of the array estimations is close to the correct value. If the conditions 
are such that none of the epochs produce an acceptable estimation, the traditional approaches of 
averaging may produce more robust results. 

b) In the case where the thermal noise at the antenna is significant and averaging over several mea- 
surements is needed, the averaging can be performed on a subset of the epochs with the smallest 
power. 

c) It should be noted that this technique can be applied to any kind of array beamforming algorithm. 



III. Simulation results 

A. In FOV interference 

This section reports on ASSC algorithm performance compared to existing techniques (classic and 
MVDR) for the example discussed in The ASSC parameters are k = 3 and ^ = 1. Figure 
shows the original image that contains a few clusters of sources. Using the classic beamformer (Figure 
(|4j))), the resulting image (classic dirty image) has wide peaks around each cluster of sources, and the 
noise is high. The ASSC classic beamformer yield a much quieter image (Figure (|4b)). The MVDR 
(Figure (|4jl)) beamformer image has higher spatial resolution than the classic beamformer (as expected). 
The ASSC MVDR beamformer (Figure (|4£)) has higher resolution than the MVDR and has the advantage 
of a quiet image. 



B. Strong out-of-FOV interference 

This section presents examples of a very strong interference 10 6 times the power of the desired sources. 
Using an East- West array with 20 antennas logarithmically spaced — 200A. Measurement was done every 
minute for a 12-hour period. The ASSC parameters used were k = 5 (out of the 719 available orientations 
) and ii). — i. Figure (pi) shows the original observed image with 6 point sources. The strong interference 
is not seen in the image (since the interferer is out of the field of view). The output of the classic and 
MVDR beamformer (dirty images) are shown in Figures (|5})) and d5b) respectively. The entire image is 
smeared with the strong interferer sidelobes. The point sources are not seen. The output of the ASSC 
MVDR beamformer is shown in Figure ([3]i). The point sources are seen clearly and the strong interferer 
sidelobes do not appear in the image at all since only correlation matrices that are affected negligibly 
from the sidelobes are selected. 



IV. Summary 

In this paper we introduced the ASSC beamformer, a method to combine rotating/many array measure- 
ments for interference suppression. The performance of the ASSC (classic and MVDR) were demonstrated 
and compared to the classic and MVDR beamformer. For interference dominant cases, the ASSC beam- 
former obtains images with higher spatial resolution and interference cancelation than either the classic 
or the MVDR beamformer. 
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(a) Original image 




(b) Classic beamformer (c) ASSC classic beamformer 




(d) MVDR beamformer (e) ASSC MVDR beamformer 

Fig. 4. Example of dirty images for a few clusters of sources. All images are plotted using the same number of contours. 
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(a) Original image 
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(b) Classic beamformer 
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(c) MVDR beamformer 



(d) ASSC MVDR beamformer 



Fig. 5. Dirty images of 6 point sources and a strong interfering source outside of the field of view with an intensity 10 6 times larger than 
the sources in the observed field of view. 
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